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A drop impacting on a solid surface deforms before the liquid makes contact with the surface. We 
directly measure the time evolution of the air layer profile under the droplet using high-speed color 
interferometry, obtaining the air layer thickness before and during the wetting process. Based on 
the time evolution of the extracted profiles obtained at multiple times, we measure the velocity of 
air exiting from the gap between the liquid and the solid, and account for the wetting mechanism 
and bubble entrapment. The present work offers a tool to accurately measure the air layer profile 
and quantitatively study the impact dynamics at a short time scale before impact. 



Drop impact on solid surfaces, beside its inherent 
beauty, has been playing an increasingly important role 
in industrial processes as diverse as ink-jet printing, spray 
cooling, and spray coating. Since it was first studied in 
1876 by Worthington [1|], the phenomenon has received 
tremendous attention from researchers, yet our under- 
standing of this subject is still far from being complete 
(see review article [2]). A challenge in studying this prob- 
lem arises from widely different time and spatial scales 
of the involved effects. Another difficulty comes from de- 
termining relevant physical parameters that govern the 
impact dynamics. For example, beside apparent param- 
eters such as the surface roughness and wettability, the 
liquid viscosity, surface tension, and density, it was re- 
cently discovered that the ambient pressure is also a cru- 
cial parameter as it dictates the splash threshold after 
impact [3]. This finding and subsequent studies [3, Isj 
suggest that the air layer between an impinging droplet 
and a solid surface may have significant effects on the 
impact's outcomes. Hence, it is essential to understand 
how the drop and the surface interact through the air 
layer. 

On the theoretical side, a mechanism of splash forma- 
tion focusing on the short time scale within which the 
drop starts being deformed has been proposed [gl, • De- 
tailed analysis and simulations have been subsequently 
developed [8, 9]. On the experimental side, the dynamics 
of droplet impact at the earliest time scale have also been 
studied; one of the most remarkable phenomena is the 
detection of entrapped bubbles under an impacting drop 
[l0-13]. The existence of these bubbles indicates that 
the drop's bottom surface is deformed before it makes 
contact with the surface. There is, however, a lack of 
detailed measurements of the air layer thickness at the 
earliest time of impact, as well as the formation of the 
entrapped bubbles. 

Here we report the first direct measurement of the 
evolution of the air layer profile between an impinging 
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droplet and a solid surface using high-speed color inter- 
ferometry. We focus on the earliest time of impact when 
the liquid has not touched the surface but starts being 
deformed due to the pressure increase in the air layer be- 
tween the liquid and the solid surface. We measure the 
air flow between the droplet and the solid surface, and 
investigate the mechanism of bubble entrapment. 
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FIG. 1. (a) Schematic of the experimental setup (not drawn 
to scale) used to study droplet impact on smooth surfaces. 
A water droplet of initial diameter Do ~ 2 mm falls on a 
glass slide of average roughness 10 nm. The bottom view is 
captured by a high-speed color camera (SA2, Photron Inc.). 
The camera is connected to a long working-distance micro- 
scope and a 5x objective to obtain a 2mm field of view, (b) 
Schematic of the air film between the drop and the glass slide 
(not drawn to scale). Light is supplied from the bottom for 
illumination; reflection of light from the upper surface of the 
glass slide and from the bottom surface of the drop causes 
interference fringes captured with the color camera, (c) An 
example of an interference pattern. 

In Fig. [1^, we show a schematic of the experimental 
setup for the present work. We generate liquid drops 
by using a syringe pump to push liquid out of a fine 
needle. The drop detaches as soon as its weight over- 
comes the surface tension and then falls on a microscope 
glass slide (Menzel microscope slide, average roughness 
^ lOnm). In our experiments, the working liquid is 
milli-Q water (density = 998kg/m^, surface tension 
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FIG. 2. (a) Reference thickness of the air film between a 
lens and a glass slide. The lens has radius of the surface in 
adjacent with the glass R — 200 mm. Inset: schematic of the 
setup used to calibrate colors, (b) Color variation in the radial 
direction of the interference pattern used to calibrate colors, 
(c) Relation between thickness and reference colors. 

cr^ = 72 X IQ-^N/m, viscosity = IQ-^mVs). The 
drop typically has diameter D ^ 2 mm and its velocity 
before impacting the surface can be adjusted by vary- 
ing the needle's height H. We capture the drop impact 
from the bottom with a color high-speed camera (SA2, 
Photron Inc.) connected to a long working-distance mi- 
croscope (Navitar Inc.) and a 5X objective. The field of 
view achieved by this combination is 2 mm. We illumi- 
nate the impact area from below by supplying white light 
from a high-intensity fibre lamp (Olympus ILP-1) to the 
microscope's coaxial light port. When a drop approaches 
the glass slide, a thin film of air is formed between the 
liquid and solid surfaces before wetting occurs. Light of 
the same wavelength coming from the bottom, upon re- 
flection from both surfaces of the film (Fig. [TJd) forms 
interference patterns recorded by the camera. Each one 
of these patterns consists of constructive (bright) and de- 
structive (dark) fringes; the fringe spacing depends on the 
air layer thickness and the wavelength of incident light. 
Since the lamp emits light of multiple wavelengths, the 
superposition of all available patterns produces concen- 
tric rings of rainbow colors as shown in Fig. [It. 

In most of our experiments, we set the camera's frame 
rate to 10000 frame per second (fps), and its resolution to 
512 X 512 pixels to capture droplets with impact velocity 
less than 0.5 m/s. In the case that the impact velocity is 
higher, the frame rate can be set as high as 86400 fps at 
resolution 32 x 256 pixels to capture the impact dynamics. 

In order to extract the absolute thickness of the air 
layer between an impinging drop and a glass surface, we 
construct a set of reference colors that can be related to 
absolute thickness. We put a convex lens on top of the 
glass slide (see inset in Fig. [JJd) and observe the interfer- 
ence rings caused by the air film between two surfaces. 
Since the pattern consists of concentric rings of differ- 
ent colors, and the air film thickness is known at each 
radial location, each color along a line passing through 
the center of these rings is associated with a thickness 
value. In Fig. [JJd, we show the air thickness profile be- 
tween the lens and the glass slide. The color variation 
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FIG. 3. (a) Color sample of interference pattern taken at 
t — 1.67 ms after the bottom- view camera detected the drop. 
(6) Color difference in grayscale computed by Eq. [T] with 
candidate profiles shown in white solid lines, (c) Average 
color difference of candidate profiles shown in (6). {d) Solid 
lines: profiles computed from three different experiments with 
V — 0.22 m/s. Dashed line: V — 1.1 m/s. Note the extremely 
different length scales at the x— axis (mm) and /ir— axis (/xm) 
in {b) and {d). 

due to change in air thickness is obtained by taking a 
thin radial strip of 100 x 2200 pixels from an image of 
an interference pattern and then averaging colors in the 
transverse direction to reduce noise. The resulting strip 
(Fig. [2^), which has no color variation in the transverse 
direction, contains N = 2200 pixels in the x— direction 
and hence N reference colors that can be used for calibra- 
tion. Since the camera uses the sRGB model to represent 
colors, the color of each pixel i is represented by a color 
vector (i?i, Gi^Bi). The pixel's coordinate is x^, which is 
related to a value of thickness h\.. Thus, we have a set of 
reference colors (i?i, G^, Bi) for 1 < i < TV, each of which 
is associated with a reference thickness h\.. The reference 
thickness range is < /i* < 4 /im. The thickness-color 
relation is shown in Fig. [2t. 

The sRGB model, however, is generally not preferred 
when comparing colors between experiments because it 
does not decouple light intensity and color information, 
which poses a problem due to variations in illumination 
conditions such as light intensity, incident and observing 
angles. Instead, we use the CIE 1976 color model (also 
called CIELAB), a model that is most effective in decou- 
pling light intensity ([l3) (see the Supplemental material 
for details). A color in CIELAB model has three compo- 
nents: L for lightness information, and a and h for color 
information. Thus, we can separate light intensity from 
our analysis by omitting the component L. Each refer- 
ence color i after intensity decouphng is represented by 
a two-component vector (a^, 6* ) and is associated with a 
value of reference thickness h], ioi 1 < i < N and h], is 
in the range /im <h\.<A jam. 

In Fig. [3^, we show a color sample, which was taken 
along a diameter of an interference pattern under a drop 
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FIG. 4. Snapshots of interference patterns obtained during drop impact and their corresponding calculated profiles {V 
0.22 m/s, D — 2 mm). See the supplementary movie. 
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FIG. 5. (a) Dimple velocity Vdim vs. time, (b) The average 
velocity of air Vair vs. time at different radial locations Rc = 
0.1mm (upward triangles), Rc = 0.3mm (squares), Rc = 
0.6 mm (downward tringles) for an experiment with impact 
velocity V = 0.22 m/s. 



with impacting velocity V = 0.22 m/s. After converting 
the color of each pixel to CIELAB color space and de- 
coupling light intensity, we calculated the color difference 
dE'^^ for each color (a^, 6^) in the sample (1 < j < M = 
596) and each reference color (a* , 6* ) using the Euclidean 



distance: 



1/2 
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for 1 < i < N and 1 < j < M. Since each color (a^, bl) is 
associated with a coordinate xj and, recall that each ref- 
erence color (a^, 6* ) is associated with a value of reference 
thickness /i* , the color difference dE'^^ can be thought of 
as a function of and Xj. In figure [3)3, we show a plot of 
dE'^ in grayscale for 1 < i < TV and for l<j<M. The 
range of the index i translates to the range of reference 
thickness diS < hr < 4 /im, and j to —1 < x < 1mm. 
In the plot, black means dE = and hence zero color 
difference, whereas white means the largest color differ- 
ence. A vertical line at a particular value of x has all 
possible values of the film thickness at that point; the 
correct thickness value corresponds to the darkest point. 
In the case that there are multiple dark points on the 
same vertical line with insignificant difference between 
them, thickness determination is not trivial. We note 
that, however, the film profile is continuous and smooth. 
Evidently, there are only a few continuous dark lines that 
can be distinguished without any abrupt change in slope. 
In figure [3)3, we show the candidate profiles in white solid 
lines (labeled from (i) to (v)). The film thickness profile 
can be identified by considering the average color differ- 
ence (dE) along each candidate profile L: 
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where the sum is taken for ah the pixels along the profile 
L and then divide by the number of pixels {Nl). In fig- 
ure [3fc, we show {dE)^ for all profiles. The smallest color 
difference is along profile {iv) for which {dEY^ = 0.06, 
whereas the second smallest one is along profile {iii) for 
which {dEy^^ = 0.09. As a result, we conclude that pro- 
file (iv) is the air layer profile. A test case of an air 
film with a known thickness profile shows that the accu- 
racy of our method is within 40 nm (see the Supplemental 
material). To check the reproducibility we repeated the 
experiment several times and extracted the air thickness 
in each experiment at the same time. The computed pro- 
files are shown in Figs. [3li. Given the variations between 
experiments such as releasing time, drop size, surface 
properties, etc., the method gives remarkably consistent 
results. 

In figure m we show interference patterns obtained dur- 
ing drop impact {V = 0.22 m/s and D = 2 mm) and their 
corresponding thickness profiles of the air layer. We de- 
fine t = as the moment when the liquid completely wets 
the soHd surface. From the first pattern detected by the 
camera {t = —3.75 ms), it is readily seen that a dimple 
is already formed, which means that the camera did not 
capture the entire deformation process of the drop's lower 
surface, probably due to limited coherence length of the 
light source used in the present experiment. Subsequent 
profiles show that the dimple's height gradually reduces, 
while the liquid continues spreading in the radial direc- 
tion. At t = —0.42 ms, the liquid starts wetting the glass 
surface at one point along the rim of the air layer and 
then propagates to the other side; the wetting process 
happens faster along the rim where the air thickness is 
smallest and finally traps air bubbles (indicated by an ar- 
row at t = ms). For experiments done under the same 
conditions, although the time that wetting occurs varies, 
we observe the same bubble-trapping dynamics, that is, 
the wetting front propagates faster at the rim and finally 
enclose the air pocket underneath the drop. 

We note that there is a plateau (at x ~ ±0.5 mm) in 
the thickness profiles from t = —2.92 ms to t = —0.21 ms. 



This is due to the drop's oscillation as it falls down at a 
small distance from the surface (in this case H = 6 mm). 
When the drop detaches from the needle, capillary waves 
are generated and propagate to the other side. The sur- 
face deformation caused by these waves affects the dimple 
shape in addition to the pressure increase in the air layer 
under the drop. In the case that the drop is released 
from a larger height leaving sufficient time for viscosity 
to damp capillary waves, we do not observe the plateau 
in the thickness profiles. As shown in Fig.[3]i, the plateau 
is not present in the case of higher impact velocity case 
{V = 1.1 m/s, H = 66 mm). 

We now quantify the velocity at the center of the 
dimple Vdim (inset of Figs. [5^). As shown in Figs. [5^, 
Vdim is found to be very small (roughly two orders of 
magnitude smaller) as compared to the impact velocity 
V = 0.22 m/s, which implies that the fiuid at the bot- 
tom of the drop has decelerated before the camera starts 
capturing the interference fringes. Nonetheless, our mea- 
surements capture well the deceleration process of the 
lower surface of the drop from the detection point until it 
is brought to rest. Moreover, we estimate the horizontal 
velocity of air Vair based on the change in volume con- 
fined by a cylinder of radius Rc under the liquid surface 
(inset of Figs. [5^). In figure [SJd, we show Vair at several 
values of Rc- The data show a consistent increase of the 
air velocity at a given time as it gets closer to the rim of 
the air layer where the thickness is minimum. For higher 
impact velocity cases, the velocity of air is much higher 
due to the extremely thin air gap at the rim. 

In conclusion, we have used high-speed color interfer- 
ometry to measure the complete profile and its evolution 
of the air layer under an impacting drop for impact veloc- 
ity V = 0.22 m/s and V = 1.1 m/s. From the experimen- 
tal measurements, we account for the wetting mechanism 
which results in entrapment of bubbles after impact. We 
also experimentally quantify the velocity of air fiow be- 
tween the drop and the surface, as well as the velocity 
of the dimple before wetting occurs. Our results offer a 
benchmark for theories of drop impact. 
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